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Abstract 

Within self-consistent field theory we study the phase behavior 
of a symmetrical binary AB polymer blend confined into a thin film. 
The film surfaces interact with the monomers via short range poten- 
tials. One surface attracts the A component and the corresponding 
semi-infinite system exhibits a first order wetting transition. The 
surface interaction of the opposite surface is varied as to study the 
crossover from capillary condensation for symmetric surface fields to 
the interface localization/delocalization transition for antisymmetric 
surface fields. In the former case the phase diagram has a single criti- 
cal point close to the bulk critical point. In the latter case the phase 
diagram exhibits two critical points which correspond to the prewet- 
ting critical points of the semi-infinite system. Only below a triple 
point there is a single two phase coexistence region. The crossover be- 
tween these qualitatively different limiting behaviors occurs gradually, 
however, the critical temperature and the critical composition exhibit 
a non-monotonic dependence on the surface field. 

The dependence of the phase behavior for antisymmetric bound- 
aries is studied as a function of the film thickness and the strength of 
the surface interactions. Upon reducing the film thickness or decreas- 
ing the strength of the surface interactions we can change the order of 
the interface localization/delocalization transition from first to second. 

The role of fluctuations is explored via Monte Carlo simulations 
of a coarse grained lattice model. Close to the (prewetting) critical 
points we observe 2D Ising critical behavior. At lower temperatures 
capillary waves of the AB interface lead to a pronounced dependence 
of the effective interface potential on the lateral system size. 

* email: Marcus.Mueller@uni-mainz.de 

1 Introduction 

Confining a binary mixture gives rise to a rich interplay between wetting 
and miscibility behavior. |l], H J|, ij In a porous material or a slit-like pore 
with identical surfaces the criticalpoint of the mixture is shifted away from 
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its bulk value to lower temperatures and higher compositions of the com- 
ponent prefered by the walls. E| The phenomenon of capillary condensa- 
tion which occurs if the two confining boundaries are symmetric is well- 
known.^] In some applications (e.g., coatings or dielectrics), however, the 
surfaces of the film (e.g., a solid substrate and vacuum) interact very dif- 
ferently with the constituents of the mixture. In the opposite limit of anti- 
symmetric walls, i.e., one wall attracts the ^-component with exactly the 
opposite strength than the other wall the S-component, an interface local- 



ization/delocalization transition ||, 0, ||, ||, 10 1 occurs close to the wetting 
transition of the semi-infinite system. 

We study the dependence of the phase behavior of a polymer mixture 
confined into a thin film on the surface interactions and film thickness via 



self-consistent field calculations [O, [i~2|, 13 1 and Monte Carlo simulations. [14] 



Our paper is arranged as follows: in the next section we describe our mean 
field calculations and discuss the dependence of the phase diagram on the 
monomer-wall interactions. Then we investigate the effect of fluctuations 
via simulations of a coarse grained lattice model. The paper closes with a 
discussion of the anticipated chain length dependence and an outlook on 
experimental systems. 



2 Self— consistent field calculations 

We calculate the phase behavior of a confined AB mixture within the self- 
consistent field theory of Gaussian polymers, [pi], 16 The film comprises a 
volume Vq = Aq x L x L. Aq denotes the film thickness, while L is the 
lateral extension of the film. The density at the film surfaces decreases to 



zero in a boundary region of width A w according to [17] 



l-cos ™ 

;0<x<A w 

<bn(r) = I 1 ;A w <x<A ~A w (i\ 

V 2 w ' ;A -A w <x<A 

where $o denotes the ratio of the monomer density and the value p in the 
middle of the film. The thickness A of an equivalent film wit h co nstant 



monomer density p is A = Ao — A w . We choose A w = 0.15i? e ,[17| where 
R e is the end-to-end distance. Both surfaces interact with the monomer 
species via a short range potential H: 



4A 1 ^ e {l+cos(^)} 



; < x < A v 



H(x) = < ] A w <x <A$- A w (2) 

4A 2 R e (l+cos(^l^l N )} 

1 ^ W —L± -Av-A w <x<Av 

H > is attractive for the A monomers and repulsive for the B species. The 
normalization of the surface fields Ai and A2, which act on the monomers 
close to the left and the right surface, is chosen such that the integrated 
interaction energy between the surface and the monomers is independent of 
the width of the boundary region A^. 

A and B polymers contain N monomers and are structurally symmetric. 
The polymer conformations {r a (r)} determine the microscopic A monomer 

density $a(i") = J2a=o Jo d T <5 (r — r a (r)), where the sum runs over all 
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n A A polymers in the system and < r < 1 parameterizes the contour 
of the Gaussian polymer. A similar expression holds for $s(r). With this 
definition the semi-grandcanonical partition function takes the form: 

n +A/in A /2k B T -A/j.n B /2k B T , , 

Z ~ T = : Pa[t] Pb[t]6($ -$ a -$b 

n% n ^ n B \ J J V 

x exp (-p J d 3 r {x^A^B - #($A - $b)}) (3) 

where n = n A + n B and A/i represents the exchange potential between A 
and B polymers. The functional integral V sums over all conformations of 

the polymers with the statistical weight exp (—-^z lo ^ r (ht) J °^ a non ~ 

interacting Gaussian polymer. The second factor enforces the monomer 
density profile to comply with Eq.([l]) (incompressibility) . The Boltzmann 
factor in the partition function incorporates the thermal repulsion between 
unlike monomers, which is described by the Flory-Huggins parameter \i 
and the interactions between monomers and surfaces. 

In mean field approximation the free energy is obtained as the extremum 
of the semi-grandcanonical free energy functional: 

= + ln£- lnQ+ 1 / d 3 r x N$a$b - HN {$ A - $ B } 
nksT Vq V J 

- i J d 3 r { W A $ A + W b ^b} + E {$ - $A - $b} (4) 

with respect to its arguments W A , Wb, $b> S. Qa denotes the single 
chain partition function: 

Qa[W a ] = yj V^PM e-f> W A« T » (5) 

and similarly for Qb\ Q = ex.p(Afj 1 /2k B T)QA + exp(— A/i/2fcsT)Qs- 

The values of Wa, Wb, $ A , H which extremize the free energy func- 
tional are denoted by lower-case letters and satisfy the self-consistent set of 
equations 

w a (t) = X NMr)-H(r)N + ar) and <j> A (r) = (6) 

<&o( r ) = 0A(r) + </>s(r), and similar expressions for wb and <f> B - Substituting 
the extremal values of the densities and fields into the free energy functional 
(3) we calculate the free energy G of the different phases. At coexistence 
the two phases have equal semi-grandcanonical potential. 

To calculate the monomer density we employ the end segment distribu- 
tion which satisfies a diffusion equation. We expand the spatial dependence 
of the densities and fields in a set of orthonormal functions. fl6|, |l7 | This 
procedure results in a set of non-linear equations which are solved by a 
Newton-Raphson-like method. We use up to 80 basis functions and achieve 
a relative accuracy 10 -4 in the free energy. 

The phase diagram in a thin film (Ao = 0.9-R e ) with antisymmetric sur- 
face fields of various strengths is presented in Fig.l. For weak surface fields 
we find a second-order interface localization/delocalization transition. [13] 
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Close to the critical point of the bulk enrichment layers form gradually 
at each surface and stabilize an interface in the center of the film. The films 
remains laterally homogeneous in this "soft-mode" phase. Only below the 
second order transition the interface is localized at one wall and ^4-rich 
and £?-rich domains coexist laterally. Upon increasing the monomer-wall 
interaction strength the transition temperature shifts to lower values. 
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Fig.l: (a) Phase diagrams in a film with antisymmetric surface fields. The values 
of the surface fields AN are indicated in the key. For AN < 0.1425 we find a single 
critical point, while we find two critical points for larger surface fields, (a) dis- 
plays the phase diagram in the temperature-composition plane, while (b) presents 
the coexistence curves A/i Goex (xAf)- (For A7V = 0.15 the two critical points are 
indistinguishable on the scale of panel (b)). From |0. 



At large monomer-wall interaction the wetting transition is of first or- 
der and we observe also a first order interface localization/delocalization 
transition. [12, 13 The prewetting coexistence slightly above the first order 
wetting transition temperature gives rise to two miscibility gaps in a film. 
At each wall a thin and a thick enrichment layer of the prefered component 
coexist. These coexistence regions terminate in two critical points which are 
the analogs of the prewetting critical points at each wall. When we decrease 
the temperature the miscibility gaps open and form a triple point at which 
an ^4-rich phase, a phase where the interface is delocalized in the middle 
of the film, and a i?-rich phase coexist. Below the triple temperature an 
A-rich and a S-rich phase coexist. 

The two types of phase diagrams are separated by a tricritical interface 
localization/delocalization transition. For our model we find a tricritical 
transition at AiV = 0.1425 and Ao = 0.9i? e . When we increase the film 
width the tricritical strength of the surface interactions decreases. At tri- 
criticality, there is only a single critical point, but the binodals are charac- 
terized by an exponent /3 = 1/4 in mean field approximation. 

Fig. lb presents the coexistence curves as a function of temperature and 
exchange potential. For a second order or a tricritical transition phase 
coexistence occurs at A/i = due to the symmetry of the system with 
respect to exchanging A ^ B. For a first order transition the coexistence 
curve lays on the symmetry axis A^ = for low temperatures. When the 
temperature rises, it bifurcates at the triple point; each continuation is the 
thin film analog of the prewetting line at the corresponding surface. Since 
the coexisting phases do not possess the symmetry of the Hamiltonian, 
A/^coex 7^ 0. For short range monomer-wall interactions the prewetting-like 
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curves deviate linearly from the bulk coexistence value jl^] and terminate in 
critical points. 
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Fig. 2. (a) Phase diagram for A./V = 0.5 and various film thicknesses Ao. For 
Ao = 2.6i? e and 0.9i? e the interface localization-delocalization transition is first 
order, Ao = 0.605i? e corresponds to a tricritical transition, while the transition is 
second order for Ao = 0.5i? e . (b) Phase diagram as a function of temperature and 
chemical potential for the same parameters than in (a) . From |0| 



The interface localization/delocalization transition can also be changed from 
first to second order by decreasing the film thickness, g, |§ pL [jj be- 
cause the effective interactions between the interface ana each walTinterfere. 
The three minimum structure of the effective interface potential close to the 
triple point changes to a single minimum as the film thickness becomes com- 
parable to the range of the effective interaction between the interface and the 
wall. The film thickness dependence of the miscibility behavior is presented 
in Fig. 2. For Ao < 0.605i? e the interface localization/delocalization transi- 
tion is second order even though the wetting transition in the semi-infinite 
system is first order. 

Of course, surface interactions in experimental realizations are never 
strictly antisymmetric or symmetric and it is important which degree of 
asymmetry is permissible without loosing the qualitative features of the 
limiting cases. The phase diagram for non-symmetric boundary fields is 
discussed in Fig. 3. The right walls attracts the ^4-component of the mixture 
and the surface fields lead to a first order wetting transition in the semi- 
infinite system. The monomer-wall interactions at the opposite wall are 
tuned from attracting A (symmetric boundaries, capillary condensation) 
to attracting B (antisymmetric boundaries, interface localization/delocal- 
ization) . 

For symmetric boundaries the critical point is shifted to lower tempera- 
tures and higher concentration of the A species attracted by both surfaces. 
The coexistence value of the chemical potential A^ coex is shifted to values 
disfavoring the A component and the shift is roughly proportional to the 
inverse film widths (Kelvin equation). Above the wetting transition tem- 
perature an almost pure ^4-rich phase coexist with a 5-rich phase. In the 
latter phase, there are thick enrichment layers of A at the walls and the 
B component prevails at the center of the film. When the temperature is 
lowered towards the wetting transition temperature T wet , the thickness of 
the enrichment layers rapidly decreases. This gives rise to a convex curva- 
ture of the B-rich binodal slightly above T wet . Below T wet the coexisting 
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phases are almost pure , i.e., the enrichment layers in the B-rich. phase are 
negligible and A/j, cocx is independent of temperature. 
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Fig. 3. (a) Binodals for A = 2.6R e and AiN = 0.5. A 2 N varies as indicated in 
the key. The dashed curve shows the location of the critical points. Filled circles 
mark critical points, open circles/dashed horizontal lines denote three phase coex- 
istence for = —0.3675 and —0.5. The inset presents part of the phase bound- 
ary for antisymmetric boundaries, (b) Coexistence curves in the x_/V-A^ plane. 
A2^V varies according to the key. The "quasi-prewetting" lines for A/j, < and 
A2^V = —0.3675 and —0.5 are indistinguishable, because they are associated with 
the prewetting behaviour of the surface with interaction A\N = +0.5. FromjTll 



As we reduce the preference of the left wall for the A component the 
system becomes more symmetric and the critical point shifts to more sym- 
metrical composition and higher temperature (i.e., closer to the bulk critical 
point at xN = 2 and (j) = 1/2). If we make the left wall repulsive for the A 
component (attractive for the i3-component) the character of the transition 
gradually changes from a bulk-like unmixing transition, where the compo- 
sition of the two phases varies little spatially across the film but differs 
between the two phases, to a prewetting-like transition, where an interface 
runs parallel to the walls and the distance between the wall and the inter- 
face is the order parameter. As this change occurs the critical temperature 
(composition) passes through a maximum (minimum). Upon approaching 
the strictly antisymmetric limit the coexistence curve (at low temperatures) 
approaches the bulk value. When it intersects with the prewetting line of 
the left wall, which attracts the -B-species, a second two phase region opens 
between a thin and a thick enrichment layer at the left wall. 



3 Monte Carlo simulations 

The self-consistent field calculations presented in the previous section ne- 
glect fluctuations. In a cylindrical pore, fluctuations destroy a true phase 
transition. In a thin film they change the universality class of the critical 
points from 3D Ising-like in the bulk to 2D Ising-like in a film. Moreover, 
there is an strong interplay between the wetting behavior of the semi-infinite 
system and the phase diagram of a film with antisymmetric boundaries and 
interface fluctuations (i.e., capillary waves) might modify the predictions of 
the mean field theory. 
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Being interested in the universal features of the phase diagram we h ave 
investigated a coarse grained lattice model of a binary polymer blend[20] 
via Monte Carlo simulations. In the framework of the bond fluctuation 
model [El]] each monomer occupies the corner of a unit cell of a simple cubic 
lattice from further occupation. Monomers along the polymer are connected 
via bonding vectors of length 2, v6, 3 or 10 in units of the lattice spacing. 
Monomers interact through a square well potential which comprises the 54 
nearest sites on the lattice. A contact between like monomers lowers the 
energy of the system by an amount e (measured in units of /cgT), while a 
contact between different species increases the energy by the same amount. 
These interactions set a temperature scale and lead to a liquid-liquid phase 
separation at T = l/e c = 69.3(3) in the bulk. The parameter e is related 
to the Flory-Huggins parameter via \ = 2z c e, where z c = 2.65 denotes 
the numb er of monomers of other chains in the range of the square well 
potential. [20] 

We work at a chain length N = 32 and monomer number density p = 
1/16. The molecules end-to-end distance is R e ~ 17 in units of the lattice 
spacing. We study thin film of geometry L x L x A. Periodic boundary 
conditions are applied in the two lateral directions, while there are hard 
impenetrable walls a distance A apart. Monomers in the two layers nearest 
to the walls interact with the boundaries. An A-monomer close to the 
right wall decreases the energy by e w = 0.16, while a £>-monomer increases 
the energy by the same amount. For symmetrical walls the interactions at 
the left wall are identical to the right wall; for antisymmetrical walls B— 
monomers are attracted by the left wall and A-monomers repelled. Using 
the Young equation we have determined the wetting transition temperature 
to T we t = 14.1(7). pq] The wetting transition occurs well inside the strong 
segregation limit and is of first order. 

We use the semi-grandcanonical ensemble, i.e., we fix the temperature 
and the exchange potential Afi and monitor the composition cp and its fluc- 
tuations. In addition to Monte Carlo moves which update the conforma- 
tions of the polymers on the lattice (local hopping attempts and slithering- 
snake-like motions) we try to change the chain species A B at fixed 
conformation. This semi-grandcanonical simulation scheme is employed in 
junction with reweighting methods as to encourage the system to sample 
all compositions with roughly equal probability. 




Fig. 4. Probability distribution of 
the composition scaled to unit 
norm and variance for various 
film thickness A (R e — 17 in 
units of the lattice spacing) and 
temperatures as indicated in the 
key. The probability distribu- 
tion of the 2D tricritical univer- 
sality class is represented by cir- 



cles. From Q 



In the simulations we monitor the joint distribution of the composition, 
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energy and monomer-wall interaction; and we use finite size scaling tech- 
niques and histogram reweighting methods to locate the critical point. As 
an example Fig. 4, presents the probability distribution of the composition 
as a function of the film thickness in the vicinity of the tricritical point. 
The distribution shows a clear three peak structure. The temperature is 
adjusted such that the ratio of the heights of the central and outer peaks 
equals 1.2. This value corresponds to the ratio of the distribution of the 
tricritical universality class. The latter quantity has been measured in sim- 
ulations of the Blume-Capel model at its tricritical point. 24] To compare 
the distributions without adjustable parameter we scale them to unit norm 
and variance. For film thickness A = 14 ~ 0.82i? e the shape of the distri- 
bution matches closely the universal function and this holds also true for 
larger systems. For A < 0.82i? e we find a second order interface local- 
ization/delocalization transition, while there is a first order transition for 
A > 0.82i? e . 

The probability distribution of the composition yields information about 
the interaction g(l) between the walls and the interface. Since the wetting 
transition occurs in the strong segregation limit, the coexisting phases in the 
bulk are almost pure </>cocx ~ or 1 and "bulk" -like composition fluctuations 
can be neglected. Hence, the distance I between the wall and the interface is 
given by I = A(p and the effective interface potential g(l) can be measured in 
the Monte Carlo simulations according to g(l) = — ksT/L 2 lnP(^>), where 
P(<t>) denotes the probability distribution of the composition. 
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Fig. 5. Dependence of the effec- 
tive interface potential g{l)/kBT 
on the lateral film extension L for 
A = 48 w 2.8i? e and e = 0.069. 
The inset presents an enlarged 
view on the minimum of the lo- 
calized state. From Ol 



The result for the effective interface potential g(l) in the vicinity of the 
triple point is presented in Fig. 5. The three minima correspond to the A— 
rich phase, the phase with the delocalized interface, and the U-rich phase. 
Unlike the situation at the tricritical point the position of the minima does 
not depend on the lateral system size. However, the effective interface po- 
tential does depend on L; the minima broaden upon increasing the lateral 
system size L and the free energy of the delocalized state decreases with re- 
spect to the localized ones. The dependence of g(l) on the lateral system size 
gives rather direct evidence for a renormalization of the effective interface 
potential by interface fluctuations p5|, 26] in the framework of a microscopic 
model. The interface in the simulations is not ideally flat, but there are 
long wavelength fluctuations of the local interface position (i.e., capillary 
waves). Since the interface is not unconstraint, the interface potential im- 
parts a parallel correlation length onto the capillary waves. £y ~ \Jdg 2 /<X 2 4>. 
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This parallel correlation length is larger in the delocalized state than in the 
localized one. For the parameters of the simulation the lateral system size 
and the parallel correlation lengths are of the same order of magnitude. For 
lateral distances smaller than £y the local position fluctuates like a free in- 
terface, for lateral distance that exceed £y interface fluctuations are strongly 
suppressed. |2?J In the simulations the lateral system size L also cuts off in- 
terface fluctuations when L < £ii. p£|l Interface fluctuations reduce the free 
energy of the system and the minima corresponding to the delocalized state 
benefits more from an increase of the lateral system size. This rationalized 
qualitatively the effect observed in Fig. 5 and quantitative description shall 
be presented elsewhere. [14] The effect is important for accurately locating 
the triple temperature. 
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Fig. 6. Phase diagram in terms of composition and temperature (a) and exchange 

potential and temperature (b) for film thickness A48 « 2.8R e . From 



The phase diagrams for A ~ 2.8R e , as revealed by Monte Carlo simula- 
tions in Fig. 6, confirm the qualitative predictions of the self-consistent field 
calculations. For symmetric boundary conditions there is a small shift of 
the critical point to lower temperatures and higher concentration of the 
A species prefered by the walls. The "bulge" of the A-poor binodal is a 
consequence of the vicinity of the prewetting critical point. [E^] If the film 
thickness were larger the coexistence value A// coex would be smaller (Kelvin 
equation) and the coexistence curve would intersect the prewetting line. In 
this case there would be a two phase coexistence region also in the case of 
symmetric boundaries.^, ^] In the antisymmetric case the phase diagram 
comprises two critical points - the analogs of the prewetting critical points 
at each surface. The concomitant miscibility gaps joint to form a triple 
point ultimately above the temperature of the first order wetting transition 
temperature of the semi-infinite system. Unlike the mean field predictions, 
however, the shape of the binodals is much flatter reflecting the 2D Ising 
behavior close to the critical points. 



4 Discussion 

We have calculated the phase diagram of a symmetric polymer mixture 
confined into a thin film in mean field approximation and by Monte Carlo 
simulations. The mean field calculations reveal a rich interplay between 
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the phase behavior in confined geometry and the wetting behavior of the 
semi-infinite system and the general features of the phase behavior are 
confirmed by our Monte Carlo results. However, fluctuations result in two 
modifications of the mean field results: (i) In the vicinity of the critical 
point we observe 2D Ising critical behavior with much flatter binodals than 
the parabolic binodals of the mean field universality class, (ii) The effective 
interface potential is renormalized by capillary waves. This leads, e.g., to 
a systematic overestimation of the triple temperature by the mean field 
treatment. 

Qualitatively the interplay between the prewetting behavior and the 
phase diagram in a film with antisymmetric boundaries is not specific to 
polymer blends but is rather characteristic of all binary mixtures. We ex- 
pect, however, polymer mixtures, be particularly suitable model systems for 
exploring these effects experimentally. 

If we measure the incompatibility of the species/temperature by \N and 

the length scale in units of the end-to-end distance R e = by/~N, the typical 
scale of the free energy in a volume R\ is given by V~N, where where the 
reduced chain length N = (pR^/N) 2 measures the degree of interdigitation. 

The importance of fluctuations close to the critical point can be gauged 
by the Ginzburg criterium. For binary polymer blends one finds that fluc- 
tuations are important in the range |1 — \cN/xN\ ^ Gi with Ginzburg 
number Gi ~ 1/N. The importance of interface fluctuations can be de- 
scribed by the capillary parameter |]30|, j3l| 

W = ~A = f^i^Re) 5Q - ( 7 ) 

4VTO-AB AttVN 

where 1/A denotes the characteristic length of the interface potential and a 
is the interface tension. \R e and aR^/y/NksT are functions of \N only. 
Upon increasing the degree of interdigitation N, the capillary parameter 
uj ~ 1/V~N decreases. 

In the limit of long chains (or strong interdigitation) these fluctuation 
effects are suppressed and the self-consistent field theory is believed to de- 
scribe the confined blend quantitatively appropriately. For the bond fluc- 
tuation model N = 32 corresponds to N ~ 91 and this chain length of 
our coarse grained model corresponds roughly to 100-150 repeat units of a 
real polymer. Since fluctuations do not alter the qualitative mean field sce- 
nario for the rather short chains investigated by Monte Carlo simulations, 
we expect this a fortiori for experimental systems. Furthermore the wetting 
transition in binary polymer mixtures occurs far below the critical point. 
Bulk-like composition fluctuations are only of minor importance; the sys- 
tems are well describable via effective interface Hamiltonians and the effect 
of interface fluctuations is clearly observable. Additionally, the size of the 
enrichment layers is set by R e and, hence, is much larger than in atomic 
liquids. Sophisticated profiling techniques (e.g., nuclear reaction analysis, 
neutron reflectometry or secondary ion mass spectrometry) are available to 
accurately determine the thickness of wetting layers in experiments. 

Indeed, experiments on polymer mixtures were the first to investigate 
fluctuation effects in the delocalized state[^, 33] and the wetting transition 
in polymer blends has been observed experimentally. [[34 , 35 1 Our results 



imply, for instance, that ultrathin enrichment layers wit 



i a thickness I 



R e /5 be unstable in the temperature range T wet < T < T prewet . Such 
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instabilities slightly above T we t have been observed for compressible one- 
component polymer films .p6|] 

Of course, both the self-consistent field calculations and the Monte Carlo 
simulations deal with highly idealized models. The effect of architectural 
asymmetry between the constituents of the mixture or the role of long-range 
van der Waals interactions between the monomers and the surfaces has not 
be considered. 
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